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Abstract: 

Electron  energy  loss  spectra  are  presented  for  a  number  of  adsorbate 
overlayers  which  form  (2*1)  superstructures  on  the  fee ( 110)  surfaces  of  Ni 
and  Cu.  These  include  (2*1)  structures  of  H,  N2,  CO,  and  0.  Lattice 
dynamical  slab  calculations  are  used  to  analyze  the  vibrational  properties 
observed  in  the  loss  spectra  which  result  from  symmetry  selection  rules  for 
dipole  activity  of  the  surface  phonons.  The  calculated  spectral  densities 
together  with  the  symmetry  selection  rules,  are  used  to  analyze  proposed 
models  for  the  adsorbate  overlayer  structures. 
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1 .  Introduction 

The  study  of  the  vibrational  properties  of  clean  and  adsorbate  covered 
surfaces  is  an  area  of  fundamental  interest  in  surface  physics.  In  these 
studies,  the  application  of  electron  energy  loss  spectroscopy  (EELS)  has 
been  shown  to  be  a  powerful  technique.1  Information  regarding  the  surface 
structure  and  interatomic  forces  may  be  obtained  from  an  analysis  of  the 
vibrational  spectra  of  the  clean  surface,  while  the  adsorption  site 
symmetry  and  bond  distances  may  be  extracted  from  the  vibrational  modes  of 
adsorbate  covered  surfaces.  Previous  EELS  studies  of  surface  phonon  vibra¬ 
tions  on  clean  and  adsorbate  covered  surfaces  have  focused  mainly  on  the 
low  index  fcc(lOO)  and  f c c ( 111)  surfaces,  such  as  the  Ni(100),2-6  Ni(lll),7 
and  Cu(100)  surfaces. 3 »5 >8 

A  complete  description  of  the  inelastic  electron  scattering  from  these 
surfaces  requires  an  understanding  of  both  the  scattering  process  and  the 
dynamics  of  the  surface.  Lattice  dynamical  calculations  for  these  systems 
have  used  a  number  of  calcul ational  schemes  such  as  the  Green  function 
techniques,9  * 10  continued  fraction  method,11  and  finite  slabs. 12-1-+  In 
the  case  of  ordered  absorbate  overlayers,  direct  information  concerning  the 
adsorption  site  configuration  may  be  ODtained  in  some  cases,  from  a  simple 
symmetry  analysis  based  on  selection  rules  for  dipole  scattering.3*15 

In  this  paper  we  discuss  the  analysis  of  the  vibrational  spectra  of 
ordered  overlayers  on  the  fee ( 110)  surfaces  of  Ni  and  Cu,  which  were 
observed  by  inelastic  dipole  scattering.  The  fee ( 110)  surface  vibrations 
have  been  less  studied  in  comparison  to  the  other  two  low  index  surfaces. 
However,  recent  lattice  dynamical  calculations16  and  experimental  measure¬ 
ments17  show  a  rich  surface  phonon  spectrum  for  this  surface.  In  this 
paper  we  first  present  for  completeness  the  relevant  spectral  function  in 
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the  EELS  cross  section  and  its  relation  to  the  dipole  active  displacements 
of  the  surface  layers.  The  lattice  dynamical  calculations  employed  here 
use  the  finite  slab  method.  This  method  was  chosen  for  its  simplicity  and 
the  ease  in  which  changes  in  the  surface  structure  may  be  accommodated. 

The  particular  systems  analyzed  are  ordered  adsorbate  overlayers  of  H,  n2, 
CO,  and  0  which  exhibit  (2*1)  superstructures  on  the  Ni  and  Cu  (110) 
surfaces.  The  tendency  for  adsorbates  to  form  (2*1)  structures  on  the 
fee (110)  surface  results  from  the  unique  anisotropy  of  this  surface.  As 
shown  in  Fig.  1(a)  this  surface  consists  of  closed  packed  rows  separated  by 
troughs  along  the  [llQ]  direction.  In  the  case  of  oxygen  adsorption  on 
these  surfaces,  a  reconstruction  of  the  surface  takes  place. 

Results  for  the  ( 2* 1 ) H  structure  are  found  to  be  consistent  with  the 
buckled  bridge  geometry  deduced  from  previous  He  diffraction  and  LEED 
results.  For  the  ( 2 *  1 ) N 2  and  CO  structures  on  the  Ni  and  Cu  surfaces 
respectively,  the  spectral  density  calculations  are  in  agreement  with  a 
terminally  bonded  configuration.  The  spectra  observed  for  the  ( 2 *  1 ) CO 
structure  is  found  to  be  consistent  with  a  previously  proposed  tilting 
geometry  of  the  CO  molecule.  For  the  (2*1)0  structure  the  density  of 
states  are  presented  for  two  different  models  of  the  surface  reconstruct! on 
and  are  compared  to  the  measured  loss  spectra  from  the  Ni  and  Cu  (110) 
surfaces . 

2.  Experimental  Details 

The  experiments  were  performed  in  a  mul ti techni que  two  level  ultrahigh 
vacuum  system  which  is  described  in  detail  elsewhere.18  Briefly,  the  upper 
level  contains  facilities  for  sample  cleaning,  gas  dosing,  low  energy  elec¬ 
tron  diffraction  (LEED)  or  Auger  electron  spectroscopy  (AES)  and  thermal 
desorption  spectroscopy  (TDS).  The  lower  level  houses  the  electron  energy 
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loss  spectrometer,  which  consists  of  a  double-pass  127°  electrostatic 
deflector  for  the  monochromator  and  analyzer.19  The  scattering  geometry  is 
fixed  with  a  total  scattering  angle  of  120°.  The  scattering  plane 
containing  the  incident  and  scattered  electrons  is  defined  by  the  surface 
normal  and  the  [110]  crystal  direction.  All  spectra  were  recorded  in  the 
specular  direction  with  an  overall  resolution  of  5  meV  (40  cm-1) 
full-wi dth-hal f-maximum  (FWHM).  The  electron  impact  energy,  referenced  to 
the  clean  surfaces,  were  3.2  and  4.3  eV  for  Cu  and  Ni  respectively. 

The  samples  consisted  of  circular  disks  of  1  cm  diameter  and  2mm 
thick.  They  were  spotwelded  between  two  0.5  mm  diameter  wires  to  a  low 
temperature  mainipulator. 20  Ta  wire  was  used  for  the  Ni  sample  while  Pt 
wire  gave  a  secure  spotweld  for  the  Cu  sample.  Sample  cleaning  consisted 
of  Neon  ion  bombardment  (500  eV )  and  annealing  to  750  and  1050  K  for  Cu  and 
Ni  respectively.  The  sample  was  heated  resistively  and  its  temperature  was 
monitored  by  a  Chromel -A1 umel  thermocouple  which  was  spotwelded  to  the  back 
of  the  Ni  sample  and  placed  in  physical  contact  with  the  edge  of  the  Cu 
crystal . 

Research  grade  H2,  N2,  CO,  and  02  (Matheson)  were  used  for  gas 
dosing.  Ordered  (2*1)  overlayers  were  produced  by  monitoring  the  low 
energy  electron  diffraction  pattern  during  exposure  until  a  sharp  well 
ordered  (2*1)  pattern  was  observed.  The  ICED  apparatus  employed  here  is  a 
custom  made  unit  consisting  of  a  set  of  four-grid  LEED  optics  i ncorporati ng 
two  microchannel  plates  for  image  intensification.  Due  to  the  high  gain  of 
the  microchannel  plates,  adverse  electron  beam  effects  are  virtually 
eliminated  by  using  very  low  beam  current  densities  ( - 1 0 “ 1 0  A  cm_‘). 

3 .  Inelastic  Dipole  Scattering  from  Surface  Vibrations 


The  inelastic  scattering  of  low  energy  electrons  from  surfaces  can  be 
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described  by  two  scattering  mechanisms  which  can  be  conveniently  discerned 
by  the  appropriate  scattering  geometry.21  These  are,  dipole  scattering 
from  the  long  range  Coulomb  interactions  between  the  incident  electron  and 
the  electric  field  fluctuations  associated  with  the  surface  modes,  and 
impact  scattering  from  short  range  interactions  with  the  atomic  potentials 
of  the  surface  atoms.  Dipole  scattering  is  associated  with  an  intensity 
distribution  which  peaks  sharply  in  the  specular  direction,  corresponding 
to  very  small  wavevector  transfer  parallel  to  the  surface,  ($„  »  0.1  In 
contrast,  impact  scattering  is  typically  associated  with  a  broad  angular 
distribution  and  gives  a  small  contribution  in  the  specular  direction  in 
comparison  to  dipole  excited  modes.21 

The  theory  of  dipole  scattering  from  surfaces  is  well 
established.1'22'23  The  inelastic  cross  section  obtained  in  the  Born 
approximation  is  given  by  a  product  of  two  factors.  The  first  is  a 
kinematic  factor  which  depends  on  the  scattering  geometry  and  is 
independent  of  the  properties  of  the  substrate.  The  second  factor  is 
related  to  the  correlation  function  of  the  charge  density  fluctuations  in 
the  surface  region.23  In  the  application  to  vibrational  motion,  this 
factor  is  the  spectral  function,  S ( oo ) ,  for  the  dipole-dipole  correlation 
function  which  is  defined  as 

S(u>)  =  1/^1  e1a)t<uz(  t ) uz ( o ) > ,  (1) 

N  2  n 

where  uz  is  the  normal  component  of  the  total  dynamic  dipole  moment,  N  is 
the  number  of  surface  primitive  cells  and  the  bracket  denotes  a  statistical 
average . 

The  semi -i nf i ni te  crystal  may  be  regarded  as  a  system  of  planes  of 
nuclei,  with  each  plane  parallel  to  the  surface.  Let  the  position  of  each 
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atom  be  denoted  by  $(L,k),  where  L  =  (L„,LZ),  C„  is  a  vector  parallel 

to  the  surface  directed  to  a  particular  unit  cell,  Lz  denotes  a 

particular  atomic  layer  and  the  index  k  specifies  the  different  atoms 

within  the  unit  cell.  The  relation  between  the  normal  component  of  the 

total  dipole  moment  and  the  displacement  u(C,k)  of  an  atom  away  from 

equilibrium  is  given  by  the  effective  charge  field  e*(C,k)  defined  by, 

u  e*(L,k)  •  u(L,k) .  (2) 

2  L  ,k 

Due  to  the  translational  symmetry  of  the  lattice  we  have 

e*(L,k)  =  e*(lz,k).  Thus  only  displacement  fields  at  the  r-point,  ($„=  0, 

can  give  rise  to  a  non-zero  dipole  moment  u2. 

The  evaluation  of  the  spectral  function  S(w)  thus  involves  determining 
a  displacement-displacement  correlation  function.  In  the  harmonic 
approximation  the  function  S(oj)  may  be  evaluated  by  expanding  the 

displacement  field  u  in  terms  of  the  eigenstates  of  the  dynamical  matrix  D 
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The  matrix  D  is  the  Fourier  transformed  dynamical  matrix  given  by 

--  .  -  i([  k'L1  k'  ) 

D(  Q  ;L  ,k  ;  L  ,k '  )  =  )  U,K.,L . 

"  11  C„  [M( L_ ,k)M( L 


1 — _1/2  exp [ - i <5,,  • ( $( L ,k )-$( L 1 ,k ‘ ) ) ] , 
z’k  )J  (5) 


where  $  is  the  force  constant  matrix.25 

With  the  use  of  eqs.  (2)  and  (3)  the-spectral  function  SU)  can  be 

directly  related  to  a  phonon  density  of  states,  g(u, {v*( L2 ,k) } ) , 
projected  onto  a  normalized  effective  charge  field  v* ( L z , k )  as. 


Tie*2 

SM  =  — ^(1  +  nU))  g(y,  { v * { L z , k ) } ) , 
2M*w 


(6) 


where 


v*(l2,k) 


r-r*1/  - 


e 


tot 


e*  =  e*(L  ,k) 
tot  ' L£ , k  z 


e*(Lz,k) 

[M(Lz,k)]1/2 

2-1/2 


(7) 

(8) 


'tot 


;  e* ( L  , k) 2/M( l  ,k) 


Lz  ,k 


(9) 


and 


g  -j,  •  v*  ( L  ,k )  •  =  |  "  X  (L  ,k)  •  v*(  L  ,k)|2  -5(  j-.j  ).  (10) 

z  X  >  L  z  ,  k  \  z  z  l  x 

Here  n(j)  is  the  Bose-£i nstei n  distribution  factor,  e*tQtis  the 
magnitude  of  the  effective  charge  field,  and  M*  is  the  effective  mass  of 
the  dipole  active  field. 

The  evaluation  of  the  spectral  function  is  facilitated  by  considering 
the  symmetry  of  the  system.  Due  to  symmetry,  the  dynamical  matrix  factors 
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so  that  displacement  fields  belonging  to  different  symmetry  classes  are 
decoupled.  Furthermore,  the  dipole  selection  rule  states  that  only  the 
totally  symmetric  displacement  fields  at  the  f  point  of  the  SBZ  are  dipole 
active.  This  selection  rule  follows  simply  from  the  invariance  of  the 
total  dynamic  dipole  moment  for  all  operations  belonging  to  the  symmetry 
group  for  the  system.10  Thus  for  the  evaluation  of  the  spectral  function 
it  is  necessary  only  to  consider  those  displacement  fields  which  are 
totally  symmetric  at  the  f  point  of  the  SBZ.  The  H  and  CO  adsorbate 
structures  discussed  in  section  5  have  a  glide  line  symmetry  and  the  more 

general  formulation  of  the  dipole  selection  rule  is  needed:  Only  totally 

symmetric  displacement  fields  with  respect  to  the  full  symmetry  group  are 
dipole  active. 

In  the  modeling  of  the  effective  charge  field  we  make  the  assumption 
that  the  effective  charges  e*(Lz,k)  are  non-zero  only  for  the  adsorbate 
layer  and  the  outermost  substrate  layer.  This  assumption  is  based  on  the 

fact  that  the  efficient  screening  by  the  conduction  electrons  will  limit 

the  dipole  activity  to  the  outermost  surface  layers.  Additionally,  it  has 
been  shown  that  the  dipole  activity  should  be  modeled  by  the  bond  stretch 
displacements  between  the  adsorbate  and  the  substrate  atoms  when  a  detailed 
comparison  of  intensities  between  theory  and  experiment  is  required.26 

To  calculate  the  vibrational  projected  density  of  states  given  by  eq. 
(10)  we  use  the  finite  slab  method  for  the  phonon  dynamics.  This  involves 
the  diagonal ization  of  the  dynamical  matrix  in  eq.  (5)  for  a  finite  number 
of  atomic  layers.  The  correspondi ng  density  of  states  constructed  from 
eq.  (10)  by  projecting  the  eigenmodes  on  the  dipole  active  displacement 
field  will  then  consist  of  a  discrete  set  of  peaks  located  at  the 
ei genf requenc i es  .  The  limited  energy  resolution  in  measured  loss  spectra 
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is  described  by  replacing  these  peaks  by  a  normalized  Gaussian  distribution 
with  a  width  given  by  the  experimental  resolution.  For  a  typical  energy 
resolution  of  5  meV  such  a  calculated  spectra  is  found  to  be  converged  with 
slabs  having  as  little  as  25  layers.  Finer  details  in  the  phonon  spectrum 
may  be  resolved  with  thicker  slabs,  with  100  atomic  layers  easily 
accomodated  with  present  day  computers.  The  size  of  the  matrix  to  be 
diagonalized,  and  hence  the  computation  time,  can  be  significantly  reduced 
by  considering  only  those  displacement  fields  which  are  symmetry  allowed, 
according  to  the  dipole  selection  rule,  and  which  may  be  decoupled  in  the 
dynamical  matrix. 

4.  The  Bare  Surface 

A  necessary  ingredient  in  the  lattice  dynamical  calculations  for  the 
oare  surface  is  a  proper  description  of  the  interatomic  forces.  For  Cu  and 
Ni  substrates  an  accurate  description  is  achieved  by  a  nearest  neighbor 
central  force  constant  model  for  the  phonon  dynamics . 1 0 > 27  In  this  model 
there  is  only  one  parameter,  the  maximum  bulk  phonon  frequency  ,  which 
is  adjusted  to  fit  the  measured  bulk  phonon  spectra.  For  Cu  and  Ni ,  f.[) 
is  equal  to  29.7  and  36.7  meV  respectively.  More  sophisticated  models  of 
the  phonon  dynamics  for  Cu  and  Ni  exist  which  account  for  second  neighbor 
and  angle  bending  interactions." 3  Incorporation  of  these  extra  terms  in 
the  phonon  dynamics  only  affected  the  surface  mode  frequencies  a  few 
percent.  This  results  from  the  fact  that  the  angle  bending  and  second 
neighbor  force  constants  are  <  10T  of  the  nearest  neighbor  force  constant, 
thus  the  single  parameter  nearest  neighbor  central  force  constant  model  is 
adopted  in  our  calculations  for  its  accuracy  and  simplicity. 

In  modeling  the  surface  force  constants  we  use  the  truncated  crystal 
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approximation,  where  the  loss  of  coordination  of  the  surface  atoms  is  the 
only  effect  considered  in  determining  the  surface  force  constants.  Thus 
the  surface  force  constants  are  determined  uniquely  by  the  bulk  values. 

This  approximation  was  shown  to  be  adequate  in  describing  previous 
experimental  results  on  the  clean  Mi  and  Cu  (110)  surfaces.16’17  Recent 
measurements  of  surface  phonon  dispersions  on  Ni(100),  along  with  lattice 
dynamical  calculations,  have  concluded  a  20%  increase  in  the  surface 
interlayer  force  constant  with  a  1.7-3. 3%  contraction  of  the  top  layer.2 
Our  previous  results  indicate  that  such  a  large  change  in  the  surface  force 
constants  is  not  present  on  the  Mi  or  Cu  (110)  surfaces  and  that  the 
surface  force  constants  can  be  modeled  after  the  bulk  values  within  an 
accuracy  of  ±15%. 

Using  the  nearest  neighbor  central  force  constant  model,  the 
evaluation  of  the  dynamical  matrix  in  eq.  (5)  for  the  fee (110)  surface  is 
strai ghtforward.  The  force  constant  matrix  is  given  by 


f  (L  ,L  ' )  =  -  ? "  ( a  0 )  n  (L,L‘  )n  (L ,  L '  )  +  5  7  4> "  ( a  0 )  n  (L,L")n  (L,C")  (11) 

‘3  J  3  L . L '  L"  “  3 

where  in  a  description  based  on  pairwise  central  potentials  3"(a0)  is  the 

second  derivative  of  the  two  body  interaction  between  atoms  L  and  L' 

evaluated  at  their  equilibrium  position  a0  and  when  LxL1  are  nearest 

nei ghbors , 


na(L)L')  =  [$(L')  -  $(L)]a/|$(L' )  -  8(L)|.  (12) 

and  zero  otherwise.  With  nearest  neighbor  central  forces  j"(a9)  is  related 
to  the  maximum  bulk  phonon  frequency  as  i"(aQ)/M5  =  ?^/8,  where  Ms  is 
the  mass  of  a  substrate  atom. 

With  the  force  constant  matrices  defined  above  we  consider  a  finite 
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We  obtain  the  surface  matrices  in  a  similar  way  using  a  truncated 
crystal,  where  the  loss  in  coordination  of  the  surface  atoms  is  the  only 
effect  considered.  These  matrices  are  given  by 

-  -  cosrrq  0  0 

2  x 


d($)(;l,l) 


0 


0 


0 


1 


(15a) 


d($..;l,2)  =  d(3„;Lz,Lz  +  l)  , 


(15b) 


d($„;l,3)  =  d(<$„;Lz,Lz+2)  , 


2  -  cos^qx  0 


d($„;2,2) 


1 

4 


0  2 

0  0 


0 


0 

3 

2 


(15c) 


( 1 5  d ) 


The  surface  eigenmodes  and  eigenvectors  for  the  finite-slab  with  r 
atomic  layers  is  then  obtained  by  constructing  the  3r  *  3r  dynamical  matrix 
in  eq.  (5)  from  the  3  *  3  submatrices  given  in  eqs.  (14)  and  (15).  Once 
the  dynamical  matrix  is  formed,  standard  algorithms  are  used  to  solve  for 
the  eigenvalues  and  eigenvectors. 

In  considering  the  effects  of  the  substrate  vibrations  on  the 
vibrational  modes  of  adsorbate  covered  surfaces,  one  is  interested  in  the 
density  of  states  projected  on  the  surface  layer  for  the  totally  symmetric 


displacement  fields  at  those  high  symmetry  points  in  the  substrate  SBZ 
which  are  equivalent  to  f  point  in  the  SBZ.  For  the  (2*1)  superstructures 
on  the  fcc(l 10)  surface,  it  is  the  f  and  X  points  of  the  substrate  SBZ 
which  are  equivalent  to  the  f  point  of  the  overlayer  SBZ,  as  shown  in  Fig. 
1(b).  Thus  depending  on  the  adsorption  site  symmetry,  various  surface 
vibrations  at  the  f  and  X  points  may  become  dipole  active  by  coupling  to 
the  motion  of  the  adsorbate  overlayer.  The  vibrational  density  of  states 
projected  on  the  motion  of  the  outer  surface  layer  is  then  obtained  using 
the  eigenmodes  and  eigenvectors  obtained  from  diagonalizing  the  dynamical 
matrix  for  a  finite  slab  at  a  particular  value  of  ($„  as  described  above. 

The  density  of  states  projected  on  the  perpendicul ar  motion  of  the  Ni 
surface  layer  for  the  displacement  field  at  the  f  point  is  shown  in  Fig. 

2.  The  spectral  intensity  shows  a  resonance  with  a  maximum  at  23  meV  for 
the  perpendicul ar  motion  of  the  outer  surface  layer.  This  particular 
resonance  at  the  T  point  is  dipole  active  in  the  absence  of  an  absorbate 
overlayer  and  has  been  observed  in  our  previous  measurements  on  the  clean 
surface.*'  Previous  surface  lattice  dynamical  calculations  show  that  the 
resonance  results  from  a  split  off  state  in  a  pseudo  band  gap  resulting 
from  the  nonmonotonic  longitudinal  bulk  phonon  dispersion  along  the  [110] 
direction.  The  results  obtained  from  the  finite  slab  calculation,  shown  in 
Fig.  2,  are  in  agreement  with  results  obtained  from  more  elaborate  Greens 
function  calculations.*0 

The  projected  density  of  states  at  the  X  point  in  the  substrate  SBZ 
are  shown  in  Fig.  3  for  the  different  displacement  fields  which  are 
displayed.  Along  the  f-X  direction  the  modes  split  into  two  classes 
resulting  from  the  mirror  plane  symmetry.  The  odd  modes  are  polarized  in 
the  y  direction  and  are  decoupled  in  the  dynamical  matrix  as  seen  in  eqs . 
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(14)  and  (15).  Similarly  the  modes  with  x-z  polarization  are  decoupled  and 
can  be  considered  separately  in  the  lattice  dynamical  calculations.  The 
sharp  peaks  observed  at  15.5,  18.4,  and  29.8  meV  correspond  to  the  S x ( X ) , 

S 2 ( X ) ,  and  S 7 ( X )  surface  phonons  observed  in  earlier  surface  lattice 
dynamical  calculations. 29  Also  shown  in  the  spectral  densities  for  the 
surface  layer  are  the  contributions  from  the  bulk  subbands.  The  S r ( X ) 
surface  phonon  is  the  low  lying  Rayleigh  phonon,  split  off  from  the  lower 
transverse  bulk  subband,  and  is  polarized  in  the  z  direction.  S2 ( X )  is 
split  off  from  the  upper  transverse  bulk  subband  and  is  polarized  in  the  y 
direction.  Similarly  S7(X)  is  split  off  from  the  longitudinal  bulk  subband 
and  is  polarized  in  the  x  direction.  The  general  behavior  for  surface 
modes  to  be  split  off  from  bulk  subbands  results  from  the  lower  restoring 
forces  for  the  surface  atoms  due  to  the  loss  in  coordination. 

5.  Adsorbate  Overlayers  With  (2*1)  Superstructures 

The  degree  of  influence  of  the  adsorbate  on  the  spectral  densities  for 
the  surface  vibrations  depends  both  on  the  adsorption  site  and  the 
perturbation  in  the  elements  of  the  dynamical  matrix  describing  the 
substrate.  In  analyzing  the  various  (2*1)  superstructures  the  displacement 
fields  at  the  r  and  X  points,  shown  in  Figs.  2  and  3,  are  analyzed 
according  to  the  adsorption  site  symmetry.  The  effect  of  this  symmetry  on 
the  vibrational  density  of  states  is  illustrated  in  particular  for  the 
( 2 x 1 ) N 2  overlayer  on  Ni(110).  In  the  case  of  oxygen  adsorption  on  these 
surfaces,  the  density  of  states  is  calculated  for  two  different  models  of 
the  surface  reconstruction. 

(a)  (2*1)  H-Ni (110) 

The  interaction  of  hydrogen  with  Ni(110)  has  been  a  subject  of  many 
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recent  investigations  due  to  the  variety  of  ordered  structures  which  form, 
depending  on  the  H  coverage  and  sample  temperature.  At  low  temperatures 
below  "  180  K,  LEED30  and  He  diffraction31  measurements  have  shown  that 
for  coverages  below  one  monolayer  (ML)  a  sequence  of  H  lattice  gas 
structures  exist.  In  particular,  a  (2*1)  ordered  hydrogen  phase  occurs  at 
9=1  ML.  Results  from  He  diffraction  measurements,  which  are  consistent 
with  LEED  observations  suggest  that  the  hydrogen  atoms  are  located  in  a 
buckled  bridge  geometry  with  the  H  atoms  shifted  off  the  short  bridge  site 
toward  the  (111)  microfacets  as  shown  in  Fig.  4(a).  The  unit  cell  for  this 
structure  contains  two  H  atoms  and  possesses  glide  lines  parallel  to  the 
[llO]  direction.  The  glide  line  symmetry  results  in  the  absence  of 
(n/2,0),  with  n  odd,  diffraction  spots  in  th  LEED  pattern  observed  at 
normal  incidence. 

Figure  4(b)  shows  the  vibrational  spectrum  of  the  (2*1)H  surface.  Two 
losses  are  observed  above  the  maximum  bulk  phonon  frequency  at  79  and  130 
meV  which  are  in  agreement  with  previous  EELS  measurements.32  Based  on 
recent  cal cul ati ons j  5  these  vibrational  modes  may  be  ascribed  to  the 
symmetric  in-plane  bending  and  stretching  modes  respectively  of  the  H  atoms 
in  the  pseudo  short  bridge  sites  as  shown  in  Fig.  4(a).34  The  third  loss 
eature  at  24  meV  lies  below  the  maximum  bulk  phonon  frequency  which 
suggests  that  it  is  related  to  a  surface  mode  at  either  the  f  or  X  points 
in  the  substrate  5BZ. 

A  comparison  with  the  spectral  densities  in  Figs.  2  and  3  for  the  f 
and  7  points,  supports  the  i nterp reta ti on  of  the  observed  loss  feature  at 
24  meV  m  terms  of  the  excitation  of  the  surface  resonance  mode  at  the  f 
point  in  the  substrate  SBZ.  The  neasured  loss  intensity  of  this  surface 


resonance  mode  agrees  well  with  the  observed  loss  intensity  of  this  mode 
for  the  bare  surface.  In  this  situation  the  dipole  activity  of  this  mode 
is  not  introduced  by  the  presence  of  the  adsorbed  H  overlayer. 

A  consideration  of  the  symmetry  properties  of  the  surface  phonons  at 
the  i  point  indicates  that  all  three  modes,  SjtX),  S2(X),  and  S7 { )? )  are 
symmetry  forbidden  for  the  adsorption  site  configuration  in  Fig.  4(a). 

S i ( X )  and  S7(X),  for  example,  are  odd  under  the  glide  symmetry  operation 
along  the  [llO]  direction,  whereas  S2(X)  is  odd  under  reflection  symmetry 
through  a  plane  parallel  the  ( 1 10 )  plane.  Only  the  symmetric  vibration  at 
the  f  point  shown  in  Fig.  2  is  symmetry  allowed  to  be  dipole  active.  Thus 
the  observation  of  only  one  surface  mode  in  the  vibrational  spectrum  is 
consistent  with  the  adsorption  site  configuration  in  Fig.  4(a)  as  deduced 
from  He  diffraction  and  LEED  measurements . 

(b)  (2*1 )N2-Ni ( 110 ) 

The  absorption  of  N2  on  Mi ( 110)  represents  an  interesting 
chemisorption  system  where  both  commensurate  and  incommensurate  overlayer 
structures  are  formed.53  Unlike  CO  adsorption  on  this  surface,  however,  N2 
on  Ni(110)  displays  only  a  single  N-N  stretching  vibration  as  a  function  of 
coverage.30  CO  on  Mi ( 110) ,  for  example,  displays  three  different  C-0 
vibrational  frequencies  correspondi ng  to  different  adsorption  sites.37  The 
fact  that  only  a  single  N-N  vibrational  frequency  is  observed  for  the 
incommensurate  phase  for  e  >  0.6  ML35  ,  where  multiple  sites  are  expected 
to  be  occupied,  raises  the  question  as  to  whether  the  N-N  stretching 
vibration  can  be  used  to  discriminate  between  different  adsorption  sites.36 

The  vibrational  spectra  for  the  (2*1)N2  surface,  correspondi  ng  to  an 
N2  coverage  of  0.5  ML,  is  shown  in  Fig.  5(c).  The  N-N  stretching  vibration 
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is  observed  at  272  meV  while  the  N-Ni  vibration  is  observed  at  43  meV. 

These  values  are  in  agreement  with  previous  EELS  3b,38  and  IR39 
measurements.  Based  on  a  comparison  of  the  N-N  stretching  vibration  with 
infrared  spectra  of  inorganic  complexes  a  linearly  bonded  N2  configuration 
was  concluded.38  The  improved  resolution  in  the  spectra  shwon  in  Fig.  5(c) 
allows  a  comparison  with  lattice  dynamical  calculations  for  various 
adsorption  sites  to  be  made  in  the  spectral  range  below  the  maximum  bulk 
phonon  frequency.  Based  on  these  calculations,  an  on-top  adsorption  site 
configuration  is  confirmed  for  the  (2*1)N2  structure. 

In  analyzing  the  (2*1  )N2  structure,  the  density  of  states  is 
calculated  for  both  the  on  top  and  short  bridge  configuration  for  the 
adsorption  site.  A  hard  sphere  model  is  used  to  obtain  the  bond  angles, 
taking  the  covalent  radii  of  the  nitrogen  and  nickel  atoms  as  0.75  and  1.25 

4  respectively.  The  nitrogen-nickel  force  constant,  <}>",  is  then  adjusted 

o 

to  give  agreement  with  the  measured  N-Ni  vibrational  frequency  at  43  meV. 

For  the  short  bridge  site  configuration  the  substrate  modes  which  are 
symmetry  allowed  to  be  dipole  active  are  the  displacement  fields  which 
correspond  to  the  S7(X)  surface  phonon  and  the  surface  resonance  mode  at 
r.  The  density  of  states  is  obtained  as  described  previously  by 
constructing  the  dynamical  matrix  for  a  finite  slab.  The  unit  cell  now 
contains  2  atoms  per  cell  for  the  substrate  due  to  the  larger  unit  cell 
defined  by  the  nitrogen  overlayer.  The  submatrices  in  eq.  (14)  are  then 
6*6  instead  of  3  *  3.  Due  to  the  point  group  symmetry  of  the  short 
bridge  site,  the  displacement  fields  with  y  polarization  are  decoupled  in 
the  dynamical  matrix.  This  reduces  the  size  of  the  submatrices  which  need 
to  be  considered  to  4  *  4.  The  resulting  matrices  for  the  surface  layer  at 
5.i  =  0  are, 
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where  d( 0 ,0 )  represents  the  adsorbate  layer  and  Ms  and  MA  are  the 
substrate  and  adsorbate  masses  respectively.  0  is  the  bond  angle  between 
the  line  joining  the  two  bonding  substrate  atoms  the  line  joining  the 
substrate  and  absorbate  atoms,  ad  is  the  perturbation  of  the  dynamical 
matrix,  ad  =  d-dj,  where  d0  is  the  unperturbed  dynamical  matrix  of  the 


substrate . 
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A  detailed  comparison  of  the  relative  intensities  in  the  density  of 
states  with  the  measured  EELS  spectra  is  dependent  on  the  projection  for 
the  dipole  active  field.  Previous  analysis  of  EELS  data  for  some  adsorbate 
structures  suggest  that  the  dipole  activity  is  mainly  due  to  bond  stretch 
displacements  between  the  adsorbate  and  the  substrate  atoms.26  For  the 
short  bridge  site  the  normalized  effective  charge  field  can  be  expressed  in 
terms  of  the  angle  5  as 
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where  5  is  the  angle  between  the  dipole  active  displacement  of  the 
substrate  atom  and  the  surface  normal. 

The  vibrational  density  of  states  is  shown  in  Fig.  5(a)  for  the  two 

extreme  cases;  5  =  0°  ( — )  and  6  =  36.9°  ( - ).  5  =  0°  corresponds  to 

displacements  normal  to  the  surface  whereas  6  =  36.9°  corresponds  to  bond 
streten  displacements.  As  observed  in  previous  calculations26  the  normal 
displacements  tend  to  assign  too  large  an  intensity  in  the  bulk  phonon 
band,  and  a  value  of  5  closer  to  bond  stretch  displacements  is  more 
likely.  The  lower  intensity  of  modes  within  the  energy  range  of  the  bulk 
phonon  band,  when  the  dipole  active  displacements  corresponds  to  bond 
stretch  displacements  rather  than  for  surface  normal  displacements,  can  be 
understood  from  the  following  argument.  For  the  modes  in  the  energy  range 
of  the  bulk  phonon  band  the  motion  of  the  high  energy  adsorbate  mode 
follows  almost  instantaneously  the  motion  of  the  substrate  atoms  in  such  a 
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way  that  the  bond  distance  is  almost  fixed.  This  results  from  the  fact 
that  the  mass  weighted  displacement  field  corresponding  to  bond  stretch 
displacements  have  almost  the  dominating  dynamical  matrix  element.  Thus  a 
relative  projection  along  the  bond  stretch  displacements  gives  a  smaller 
contribution  from  the  surface  vibrations.  This  effect  is  increasingly  more 
pronounced  for  modes  with  lower  energy. 

As  predicted  by  the  symmetry  selection  rules,-  S7()0  is  observed  to  be 
dipole  active  with  an  appreciable  intensity  for  the  short  bridge  site 
configuration.  It's  frequency  is  slightly  shifted  from  29.8  meV  for  the 
bare  surface  to  27.8  meV  for  the  nitrogen  covered  surface.  The  f  resonance 
is  not  resolved  from  the  S7()0  surface  phonon  due  to  the  finite  broadening 
used  in  the  calculation.  A  comparison  of  the  experimental  loss  spectrum 
with  the  calculated  density  of  states  can  then  be  used  to  rule  out  the 
short  bridge  configuration  for  the  adsorption  site  since  there  is  no 
features  observed  at  2d  meV  in  the  measured  loss  spectrum  where  S7(X)  is 
expected. 

For  the  on-top  site  configuration,  the  symmetry  selection  rules 
indicate  that  the  displacements  fields  corresponding  to  the  S x ( X )  surface 
phonon  and  the  surface  resonance  mode  at  F  should  be  dipole  active.  Based 
on  a  symmetry  analysis  alone,  an  incorrect  conclusion  could  be  deduced  by 
analyzing  the  measured  spectra  where  no  sharp  features  are  observed  below 
the  maximum  bulk  phonon  frequency.  A  lattice  dynamical  calculation  is  then 
necessary  to  obtain  both  the  observed  frequencies  and  intensities  of  the 
surface  modes  which  should  show  dipole  activity. 

For  the  on-top  configuration  only  z  displacements  need  to  be 
considered  in  the  dynamical  matrix  since  both  x  and  y  displacements  are 
decoupled.  The  surface  matrices  for  the  on-top  site  are  then  given  by, 
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d( 0 ,0)  = 


d(0,l)  = 
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(18b) 


Ad( 1,1)  = 


(13c) 


The  density  of  states  for  the  on- top  configuration  is  projected  on  the 
relative  normal  displacements  of  the  nitrogen  and  bonding  substrate  atom  as 
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and  is  shown  in  Fig.  5(b).  From  the  calculated  density  of  states  it  is 
observed  that  both  the  Sp(^)  and  the  f  surface  resonance  are  not  observed 
with  any  appreciable  intensity  even  though  they  are  symmetry  allowed  to  be 
dipole  active.  A  small  contribution  from  bul<  states  is  obse^vl  around  :! 
meV.  The  small  contribution  from  the  $•(?)  surface  phonon  i"d  :*e 
surface  resonance  results  from  the  fact  that  the  dipole  pro  >>  .t  ■  t"o 

terminally  bonded  N2  is  along  ".he  bond  stretch  displacement.  , 

lower  intensity  for  the  surface  vibrations  for  the  reasons  ir„.  ,  sod 
previously  for  the  short  bridge  adsorption  site.  The  lack  of  any  s"i'-p 


features  within  the  bulk  phonon  bond  in  the  measured  vibrational  spa-  trp 
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for  the  ( 2 x 1 ) N 2  overlayer  is  thus  consistent  with  an  on-top  configuration 
as  seen  in  the  calculated  density  of  states.  This  conclusion  is  in 
agreement  with  the  configuration  as  deduced  empirically  from  the  value  of 
the  N-N  stretching  vibration. 

(c)  (2*1 ) CO -C u ( 110) 

The  adsorption  of  CO  on  Cu ( 110)  possesses  many  similarities  with  the 
adsorption  of  N2  on  the  Ni(110)  surface.  These  similarities  have  been 
noticed  in  previous  studies35*36  where  CO  also  forms  a  (2*1)  ordered 
phase  at  9  =  0.5  ML  and  a  compressed  i ncormensurate  phase  at  higher 
coverages/3  Similarly,  there  is  only  one  main  CO  vibrational  band  at  2088 
cm**  which  shifts  *  6  cm**  with  increased  coverage  as  observed  by  IR 
spectroscopy/*  A  second  weaker  CO  vibrational  mode  was  observed  at  2104 
cm**  at  low  coverages  and  has  been  attributed  to  CO  adsorbed  in  defect 
si  tes  /2 

The  similarities  between  the  H2/Ji(110)  and  CO/Cu ( 110)  systems  allow 
the  previous  results  for  the  (2*1)N2  system  to  be  applied  to  the 
(2*1 ) CO/Cu ( 110)  by  scaling  the  phonon  frequencies  occurring  within  the  bulk 
phonon  bands  by  0.8,  which  is  simply  the  ratio  of  the  maximum  bulk  phonon 
frequencies  for  Cu  and  Hi.  The  measured  vibrational  spectrum  for  the 
( 2 <  1 ) C 0  overlayer  is  shown  in  Fig.  6.  The  C-0  stretching  vibration  is 
observed  at  259  meV  (2089  cm**)  in  agreement  with  previous  I R* 0  *  *  1  and 
CELS'*3  measurements.  Also  observed  is  the  C-Cu  stretching  vibration  at  43 
meV.  We  did  not  observe  a  band  at  85  meV  which  was  reported  in  the 
previous  EELS  study1*3  and  is  most  likely  due  to  some  impurities  since  it 
was  only  observed  at  low  exposures.  As  seen  in  Fig.  6  no  spectral  features 

are  observed  within  the  bulk  phjnon  bands,  similar  to  the  ( 2 *  1 ) N 2  case. 
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2y  scaling  the  Ni  results,  the  S?(^)  surface  phonon  should  be  observed  at 
-  22  meV  for  CO  occupying  a  short  bridge  site  configuration  on  Cu.  Due 
to  the  null  result  and  the  value  of  the  C-0  stretching  vibration,  an  on-top 
adsorption  site  can  be  concluded  for  the  {2*1)  overlayer  structure. 

i  d)  i 2*1) CO-Ni (110) 


The  interaction  of  CO  with  Ni(110)  has  been  a  subject  of  recent 
interest  due  to  the  different  structural  conf igin  ations  which  form  on  this 
s-.rface  depending  on  the  CO  coverage.37  Previous  EELS  measurements  have 
’certified  three  C-0  stretching  vibrations  corresponding  to  different 
a  :s.;rpti  or.  sites."*'*''  Below  saturation  coverage  the  two  C-0  modes, 

«"i;n  j c c j r  at  242  ana  256  meV,  have  been  assigned  to  short  bridge  and 
c^-top  adsorption  site  configurations  respectively.  At  saturation  coverage 
corresponding  to  ■:  =  i  ML,  CD  forms  an  ordered  (2*1)  overlayer  structure, 
'he  lEEj  pattern  contains  missing  (n/2,0),  with  n  odd,  diffraction  beams, 
*nich  indicates  a  glide  line  symmetry  along  the  [llO]  direction  similar  to 
the  previously  described  hydrogen  structure  .■* 5 

The  vibrational  spectrum  for  the  ( 2 « 1 ) CO  overlayer  is  shown  in  Fig. 
7(b).  hhe  C-0  and  C-Ni  stretching  vibrations  are  observed  at  245  and  52 
neV  respectively,  in  agreement  with  previous  EELS  results.*'*  Also  observed 
in  Fig.  7(b)  is  a  loss  feature  within  the  bulk  phonon  bands  at  30  meV. 

Cased  on  LEED  and  EELS  measurements  an  adsorption  site  configuration 
consisting  of  al ternati vely  tilted  CO  molecules  in  a  short  bridge  site  has 
been  proposed'  as  shown  in  Fig.  7(a).  Recent  electron  stimulated 
desorption  <on  angular  distribution  (ESOIAD)  measurements  have  confirmed 
the  tilting  of  the  CO  molecule  in  the  (2*1)  structure  with  a  tilt  angle  of 
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17  -  20°. A  similar  tilted  geometry  has  been  determined  for  the 
(2*1 ) CO/Pt ( 110)  structure,  although  with  an  on-top  adsorption  site 
consistent  with  the  higher  CO  stretching  frequency  on  that  surface."17 

The  symmetry  analysis  pressented  for  the  (2*1 ) H  system  can  be  applied 
to  the  adsorption  site  configuration  in  Fig.  7(a)  due  to  the  similar  point 
group  symmetry.  As  was  found  for  the  hydrogen  case,  displacement  fields 
correspond!' ng  to  the  surface  modes  at  the  X  point  shown  in  Fig.  3  are 
symmetry  forbidden  to  be  dipole  active,  whereas  the  resonance  mode  at  the  f 
point  is  symmetry  allowed.  Furthermore,  two  vibrational  modes 
correspondi ng  to  the  frustrated  rotational  and  translational  modes  of  the 
molecule  will  be  dipole  active  due  to  the  lower  symmetry  of  the  tilted 
molecule.  Based  on  cluster  cal  cul  ati  ons ,  "t3  the  frequency  of  the  frustrated 
translation  is  expected  to  be  very  low,  <  10  meV,  since  to  first  order  this 
mode  does  not  involve  any  straining  of  the  nearest  neighbor  bonds.  The 
frustrated  rotation  should  be  in  the  frequency  range  of  the  metal-carbon 
vi brati on . 

Lattice  dynamical  calculations  for  1  ML  of  untilted  CO  molecules  on 
the  Ni(110)  surface  show  a  negligible  intensity  of  the  ?  resonance,  similar 
to  the  effect  observed  for  the  terminally  bonded  (2*1)N2  in  Fig.  5(b). 
rurther  tne  inherent  dipole  activity  of  this  surface  resonance  is  much 
weaker  than  the  dipole  activity  of  the  frustrated  translation  of  N2  against 
the  surface  as  evidenced  by  Figs.  4(b)  and  5(c).  This  leaves  the 
assignment  of  the  30  meV  loss  feature  in  Fig.  7(b)  to  be  due  to  the 
frustrated  rotational  mode  of  the  tilted  CO  molecule.  Similar  dipole 
active  losses  have  been  observed  and  assigned  to  frustrated  rotations  for 
tilte  CO  molecules  on  other  transition  metal  surfaces. "t3  The  proposed 
short  bridge  site  geometry  in  Fig.  7(a)  is  also  consistent  with  the  value 
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of  the  C-0  stretching  vibration  observed  at  245  meV.  For  a  tilted  on-top 
adsorption  site  configuration,  an  additional  mode  corresponding  to  the 
S 2 ( ^ )  surface  phonon  is  symmetry  allowed  to  be  dipole  active.  The  lack  of 
any  additional  features  at  ~  18  meV  expected  for  the  S2(50  surface  phonon 
thus  lends  additional  support  for  the  short  bridge  site  configuration  shown 
in  Fig.  7(a). 

(e)  (2*1 ) 0 - N i ( 110) ,  C u ( 110) 

The  interaction  of  oxygen  with  the  nickel  and  copper  (110)  surfaces 
has  been  the  subject  of  many  recent  structural  investigations. 

Considerable  interest  has  been  attached  to  the  (2*1)  structure  which  forms 
on  both  these  surfaces  and  shows  an  oxygen  inouced  reconstruct!' on  of  the 
surface  layer.  For  the  Cu  surface,  low  energy  ion  scattering50  and  He 
diffraction3'  have  concluded  a  missing  row  model  for  the  surface 
reconstruct! on  where  every  [001]  row  is  absent.  The  Ni(110)  surface  has 
received  considerably  more  attention  with  conflicting  views  on  the 
particular  model  for  the  surface  reconstruct! on .  Two  of  the  most  favored 
models  for  the  Ni  reconstruct!' on  are  the  missing  row  (MR)  model52"55  and 
the  sawtooth  (SW)  model57’38  which  are  shown  in  Fig.  8.  The  two  models 
differ  by  an  additional  missing  row  in  the  second  layer  for  the  SW  model. 

It  would  be  expected  that  the  surface  vibrational  modes  would  be  very 
sensitive  to  such  a  loss  of  coordination  of  the  surface  atoms  and  can 
therefore  serve  as  a  means  of  distinguishing  between  different  models  of 
the  surface  reconstruct! on . 

In  the  present  paper  we  present  results  from  EELS  for  the  (2*1) 
reconstruct! on  on  both  the  Cu  and  Ni  (110)  surfaces.  The  surface  modes 
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appearing  within  the  bulk  phonon  bands  are  expected  to  be  sensitive  to  the 
particular  reconstruction  reflecting  the  additional  loss  of  coordination  of 
the  surface  atoms  in  the  reconstructed  geometry.  Lattice  dynamical 
calculations  are  shown  for  the  two  different  models  of  the  reconstruct!' on 
for  Ni(llO).  The  results  for  Cu  can  be  obtained  by  simply  scaling  the 
results  by  the  ratio  of  the  maximum  bulk  phonon  frequencies. 

The  energy  loss  spectrum  for  the  (2*1)0  surfaces  are  shown  in  Fig.  9. 
Nearly  identical  oxygen-metal  stretching  vibrations  are  observed  at  48  and 
49  meV  on  the  Cu  and  Ni  surfaces  respectively.  These  values  are  in 
agreement  with  previous  EELS  measurements  on  these  surfaces.55’56 
Additional  sharp  losses  are  observed  below  the  maximum  bulk  phonon 
frequency  at  15  and  30  meV  for  Ni  and  24.5  meV  for  Cu.  Both  energy  loss 
and  gain  peaks  are  observed  with  an  intensity  ratio  which  is  determined  by 
the  3ose-Einstein  distribution  factor  and  the  crystal  temperature.  In 
comparing  the  EELS  spectra  on  the  two  surfaces,  it  is  observed  that  the 
ratio  of  the  30  to  24.5  meV  loss  scales  as  the  ratio  of  the  maximum  bulk 
phonon  frequencies.  This  suggests  the  possibility  that  these  vibrations 
are  related  to  the  motion  of  the  metal  substrate  atoms  and  that  the 
reconstruct! on  is  similar  on  both  surfaces.  The  15  meV  loss  observed  on 
the  Ni  surface  would  then  be  shifted  down  to  -  12  meV  on  the  Cu  surface 
where  it  would  not  be  resolvable  from  the  elastic  beam. 

In  calculating  the  density  of  states  for  the  dipole  active  losses,  we 
use  the  structural  models  shown  in  Fig.  8.  The  position  of  the  oxygen  atom 
has  been  determined  for  the  Ni(llQ)  surface  by  ion  scatteri  ng5'*  < and 
helium  di f fracti on5 3  measurements  to  be  in  the  long  bridge  site  0.25  A 
above  the  center  of  the  top  most  nickel  rows.  With  this  geometry  the  0- Ni 
bond  distances  are  1.77  and  1.95  A  for  the  first  and  second  layer  Ni  atoms 
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respectively.  The  0-Mi  force  constant  with  the  second  layer  Mi  atoms  would 
be  expected  to  be  comparable  to  that  for  the  first  layer  Mi  atoms  due  to 

the  similar  bond  distances.  To  minimize  the  number  of  parameters  in  the 

calculation  the  0-Mi  force  constants  are  taken  to  be  the  same  for  the  first 

and  second  layer  Mi  atoms;  it's  value  is  then  uniquely  determined  by 

fitting  to  the  measured  O-Ni  vibrational  frequency  at  49  meV.  Similar  to 
the  previous  calculations,  the  Mi-Mi  force  constants  in  the  surface  region 
are  taken  to  be  equal  to  their  bulk  value. 

The  calculated  density  of  states  are  shown  in  Fig.  10(a)  for  the  MR 
model  of  the  bare  Mi  surface.  Included  in  the  results  for  the  bare  surface 
is  the  sum  of  the  density  of  states  for  the  totally  symmetric  displacements 
with  respect  to  the  long  bridge  site  for  the  first  and  second  Mi  layers. 

The  two  higher  energy  peaks  at  21.3  and  31.7  meV  result  from  a  breathing 
motion  in  the  second  layer  with  a  polarization  in  the  x  direction.  The 
lower  energy  peak  at  16.9  meV  results  from  the  perpendicul ar  motion  of  the 
top  layer  Mi  atoms. 

Shown  in  Fig.  10(b)  is  the  density  of  states  for  the  oxygen  covered 
surface.  The  projection  is  taken  on  the  perpendicul ar  motion  of  the  oxygen 
atom  to  avoid  the  complexity  in  modeling  the  effective  charge  field  for 
bonding  to  both  the  first  and  second  layer  Mi  atoms.  As  seen  in  Fig. 

10(b),  a  sharp  feature  is  observed  at  30  meV  which  is  shifted  down  slightly 
from  31.7  meV  on  the  bare  surface.  The  remaining  sharp  structure  in  the 
density  of  states  of  the  bare  surface  is  not  seen  to  be  reproduced  on  the 
oxygen  covered  surface  except  for  a  broad  feature  at  24  meV.  In  comparing 
to  the  measured  loss  spectrum,  the  30  meV  loss  feature  is  accounted  for  in 
the  calculation  for  the  MR  Model,  however  no  structure  is  observed  at  15 
neV  in  the  calculated  density  of  states  to  correspond  to  the  observed 
feature  in  the  measured  loss  spectrum. 
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The  calculated  density  of  states  for  the  SW  model  is  shown  in  Fig. 
11(a)  for  the  bare  surface.  The  lower  symmetry  in  the  SW  model  increases 
the  number  of  symmetry  allowed  diplacement  fields.  In  comparison  to  the  MR 
model  a  substantial  amount  of  intensity  is  distributed  to  lower  energies. 
This  results  simply  from  the  fact  that  the  restoring  forces  for  the  surface 
atoms  are  smaller  in  this  model  than  for  the  MR  model  resulting  from  the 
lower  coordination  of  the  surface  atoms  in  the  SW  model  in  comparison  to 
the  MR  model.  The  results  for  the  oxygen  covered  surface  is  shown  in  Fig. 
11(b)  for  a  projection  on  the  perpendicular  motion.  Sharp  spectral 
features  are  observed  at  15.4,  25.3,  and  27.8  meV.  The  features  at  25.3 
and  27.8  meV  would  not  be  resolved  with  the  5  meV  experimental  resolution 
and  would  thus  appear  as  one  loss  peak  at  -  27  meV.  Also  shown  in  Fig. 
11(b)  is  the  projection  on  the  parallel  motion  of  the  oxygen  atom  in  the  x 
direction.  Due  to  the  loss  of  the  additional  second  layer  Mi  atom,  the 
perpendi cul ar  and  parallel  motions  are  now  coupled  as  seen  in  Fig.  11(b). 

In  comparing  to  the  measured  loss  spectra,  the  SW  model  gives  two 
losses  at  approximately  the  correct  frequencies,  which  is  in  better 
agreement  than  the  MR  model.  A  decisive  distinction,  however,  cannot  be 
made  since  the  agreement  with  the  two  models  could  be  made  to  change  by 
varying  the  parameters  in  the  calculation.  The  results  do  show  that  the 
surface  modes  are  extremely  sensitive  to  the  coordination  of  the  atoms  in 
the  models  for  the  reconstruction  and  that  a  larger  data  base  obtainable  by 
angular  distribution  measurements  throughout  the  S3Z,  coupled  with  lattice 
dynamical  calculations,  could  contribute  to  the  determination  of  the 
structure  for  the  reconstructed  surfaces. 

6.  Conclusion 


A  variety  of  adsorbates  are  observed  to  form  (2*1)  superstructures  on 
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the  f c c ( 110)  surfaces  of  Ni  and  Cu.  A  lattice  dynamical  analysis  of  the 
energy  loss  spectra  of  these  adsorbate  structures  was  achieved  using  a 
simple  finite  slab  method,  together  with  the  symmetry  selection  rules  which 
are  operative  in  dipole  scattering.  For  the  (2*1)H  structure  on  Ni ( 110) 
the  observation  of  a  single  resonance  mode  at  24  meV  was  found  to  be  in 
agreement  with  the  buckled  bridge  geometry  as  deduced  from  He  diffraction 
and  LEED  measurements.  Based  on  a  comparison  of  the  loss  spectra  with 
spectral  density  calculations,  a  terminally  bonded  configuration  was 
deduced  for  the  (2*1)N2  and  CO  structures  on  the  Ni  and  Cu  surfaces 
respectively.  For  the  (2*1)C0  structure  on  Ni(110),  a  loss  feature  was 
observed  at  30  meV  which  is  assigned  to  the  CO  frustrated  rotational  mode 
for  the  proposed  structure  consisting  of  alternating  tilted  CO  molecules  in 
the  short  bridge  site.  In  the  case  of  oxygen  adsorption  on  these  surfaces 
a  reconstruction  of  the  surfaces  takes  place.  The  vibrational  density  of 
states  of  the  reconstructed  surface  was  calculated  for  the  missing  row  and 
sawtooth  models.  For  the  sawtooth  model  a  significant  amount  of  intensity 
is  redi stri buted  to  low  energies  due  to  the  lower  coordination  of  the 
surface  atoms.  A  comparison  of  the  measured  loss  spectra  with  the  density 
of  states  for  the  oxygen  covered  surfaces  favors  the  sawtooth  model  for  the 
surface  reconstruction .  However,  a  firm  distinction  between  the  two  models 
requires  a  larger  data  base  from  measurements  throughout  the  surface 
Brillouin  zone. 
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Figure  Captions 

Fig.  1  Atomic  arrangement  of  the  fcc(llO)  surface  and  the  corresponding 
surface  Brillouin  zone  (SBZ).  The  dashed  line  shows  the  reduced 
SBZ  for  an  adsorbate  overlayer  in  a  (2x1)  configuration. 

Fig.  2  Calculated  vibrational  density  of  states  for  the  perpendi cul ar 

motion  of  the  outermost  surface  layer  at  the  f  point  for  the  clean 
Ni  surface.  The  peak  at  23  meV  corresponds  to  a  vibrational 
resonance  of  the  outermost  surface  layers.16*17  The  calculation 
used  a  130  layer  slab  with  1  meV  Gaussian  broadening. 

Fig.  3  Calculated  vibrational  density  of  states  for  the  clean  Ni  surface 
at  the  X  point  in  the  substrate  SBZ.  (a)-(c)  Projections  on  the 
notion  of  the  outermost  layer  corresponding  to  the  depicted  dis¬ 
placement  fields  at  X.  The  sharp  peaks  correspond  to  the  surface 
phonons  S^tX),  S2(X),  and  S7(X),  according  to  the  notation  in 
Ref.  29.  The  calculations  were  performed  on  a  40  layer  slab  with 
2  meV  Gaussian  broadening. 

Fig.  4  (a)  Model  for  the  (2xl)H  structure  on  Ni ( 110) .  The  smaller  filled 

circles  represent  hydrogen  atoms,  (b)  Electron  energy  loss 
spectrum  for  the  (2xl)H  overlayer  on  Nit  110)  at  80  K. 

Fig.  5  A  comparison  between  results  from  lattice  dynamical  calculations 

and  the  experimental  loss  spectrum  for  the  (2x1)  ordered  overlayer 
on  Ni ( 110) .  (a)  Spectral  density  for  the  short  bridge  adsorption 

site  configuration  for  two  different  dipole  projections;  ( - ) 

projection  on  the  bond  stretch  displacements,  and  ( — )  projection 
on  the  surface  normal  displacements,  (b)  Spectral  density  for  the 
on  top  adsorption  site  configuration,  (c)  The  experimental  loss 
spectrum  for  the  Ni ( 1 10 ) - ( 2x1 )N2  overlayer  at  80  K.  Note  the 
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break  in  thr  energy  scale  to  display  the  N-N  stretching  vibration 
at  272  meV . 

Fig.  6  Electron  energy  loss  spectra  of  Cu(ilO)  with  a  (2x1 )C0  overlayer 
at  80  K. 

Fig.  7  (a)  A  schematic  model  for  the  Ni ( 110) -( 2x1 )C0  structure.  The 

smaller  filled  carbon  and  open  circles  (oxygen),  represent  the 
tilted  CO  molecule,  (b)  Electron  energy  loss  spectra  of  Ni ( 110) 
with  a  ( 2x1 ) CO  overlayer  at  80  K. 

Fig.  8  Models  of  the  (2x1)  reconstructed  (110)  surfaces  formed  by  oxygen 
adsorption,  (a)  Missing  row  model  (b)  Saw  tooth  model.  The  snail 
filled  circles  represent  oxygen  atoms.  The  numbers  denote 
sequential  layers  of  the  substrate  with  number  one  corresponding 
to  the  surface  layer.  For  each  model  both  the  top  and  side  v'ews 
are  shown. 

Fig.  9  Electron  energy  loss  spectra  of  (2x1)  oxygen  over!  ayes  .  "  ....  inn 
Nit  110)  at  300  K. 

Fig.  10  Results  from  lattice  dynamical  calculations  for  tre  mi g  ■  * 
model  of  the  surface  reconstruction  shown  in  Fig.  •  a ) .  a 
summation  of  the  vibrational  density  of  states  fpr  projections  ;.o 
the  totally  symmetric  displacements  with  respect  to  the  ’em 
bridge  site,  for  the  first  and  second  layers  of  the  bare  M 
surface,  (b)  Vibrational  density  of  states  projected  on  the 
perpendi cul ar  motion  of  the  oxygen  atom  in  the  long  bridge  site 
for  the  missing  row  model.  The  calculations  used  a  99  layer  slab 
with  a  1.5  meV  Gaussian  broadening. 
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pig.  11  Results  from  lattice  dynamical  calculations  for  the  sawtooth  model 
of  the  surface  reconstruct! on  shown  in  Fig.  3(b).  (a)  A  summation 

of  the  vibrational  density  of  states  for  projections  on  the  total¬ 
ly  symmetric  displacements  with  respect  to  the  long  bridge  site, 
“'or  the  first  and  second  layers  of  the  bare  Ni  surface,  (b) 

tibrational  density  of  states  projected  on  the  perpendi cul ar  ( - ) 

and  parallel  (-  -  -)  motion  of  the  oxygen  atom  in  the  long  bridge 
site  of  the  sawtooth  model.  The  calculations  used  a  99  layer  slab 
«i tn  a  1.5  meV  Gaussian. 
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